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ABSTRACT 

We present a study of the intrinsic deprojected ellipticity distribution of the satellite 
dwarf galaxies of the Andromeda galaxy, assuming that their visible components have 
a prolate shape, which is a natural outcome of simulations. Different possibilities for 
the orientation of the major axis of the prolate dwarf galaxies are tested, pointing 
either as close as possible to the radial direction towards the centre of Andromeda, or 
tangential to the radial direction, or with a random angle in the plane that contains 
the major axis and the observer. We find that the mean intrinsic axis ratio is ~ 1/2, 
with small differences depending on the assumed orientation of the population. Our 
deprojections also suggest that a significant fraction of the satellites, ~ 10%, are tidally 
disrupted remnants. We find that there is no evidence of any obvious difference in the 
morphology and major axis orientation between satellites that belong to the vast thin 
plane of co-rotating galaxies around Andromeda and those that do not belong to this 
structure. 
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1 INTRODUCTION 


The Local Group can be thought of as the closest labo¬ 
ratory for testing cosmological predictions on small scales, 
including the dynamical formation and evolutionary history 
of galaxies ( Freeman fc Bland-Hawthorn|2002 I. In this con¬ 
text the dwarf spheroidal (dSph) galaxies are of particular 
interest, as they represent the smallest of galactic systems, 
and are thought to be embedded in primordial dark mat¬ 
ter subhalos, thus giving us insights on the earliest epochs 
of galaxy formation in the Universe. The investigation of 
these satellites is thus one of the most active areas of Local 
Group research at present ( Tolstoy et al.|[2C)09 |. Their mor¬ 
phological parameters, their masses, their spatial distribu¬ 
tion around their host, or the number of observed satellites, 
can be confronted with predictions of the A Cold Dark Mat¬ 
ter (ACDM) model ( [Kravtsov et al.||2004| [Tollerud et al.j 
2014). For example, their spatial distribution can provide 


insights on the shape of the dark matter halo of the host 
galaxy (Penarrubia et al. 2002), whilst the dwarf galaxies 


morphological characteristics can provide insights on their 
formation history ( Zolotov et al.|2011| Carraro|2014 ). 

The very nature of the dSphs as faint, poorly populated 
systems, poses several observational challenges, even for 
their identification in the Local Group. Generally, they have 
been discovered through the identification of small overden- 
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sities of stars over the contaminating foreground stellar pop¬ 
ulations (and background population of mis-identihed com¬ 
pact galaxies). Improvements in instrumentation and anal¬ 
ysis techniques have led to the discovery of several new Lo¬ 


cal Group dSph galaxies in recent years (Belokurov et al. 
2006[ jZucker et al.||2006[ [Walsh et al.| |2007| [McConnachie 


et al.[ 2008[ 


Martin et al. 


2009[ [2013[ ). The known popu¬ 


lation of dSphs in the Local Group has thereby increased 


considerably (McConnachie 2012), which is now beginning 
to allow us to draw statistically firm conclusions about the 
properties of the dSph population. For instance, a recent 
discovery has shown that half of the dwarf galaxy satellites 
of the Andromeda galaxy (M31) reside in a vast thin plane 
around their host (VTP), and appear to have coherent mo¬ 
tions ( Ibata et al.[2013 ). This observational dynamical align¬ 
ment is still not fully understood within the ACDM frame¬ 
work where the subhalos that might host dSphs are a priori 


expected to possess only a weak alignment |Shaya & Tully 

2013| 

Ibata et al.|2014||Pawlowski et al.|2014||Bahl & Baum- 

gardt 

2014 

1 . Finding whether the satellites that partake in 


the alignment have a particular intrinsic shape and orienta¬ 
tion that distinguishes them from the other dSphs orbiting 
around M31 may help in assessing whether the two sub¬ 
populations have experienced a different dynamical history. 


The morphological properties of satellite dwarf galaxies 
can also provide tests of theoretical frameworks of galaxy 
formation on small scales. These galaxies are usually de¬ 
scribed as a baryonic component embedded in a massive 
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dark matter halo (Walker et al. 20091, meaning that their 


have a prolate shape, like a cigar. This is indeed the natural 


stars are relatively protected from the tidal field of their 
giant host. But the formation of a plane of satellites (as ap¬ 
parently required in M31), or any past interactions could 
have produced strong tidal effects that could have strongly 
affected their baryonic distribution and make them partly 
lose their dark matter halo. Another possibility is that some 
dSphs galaxies are actually not embedded in dark matter 
halos at all, but are rather remnants of so-called tidal dwarf 
galaxies which would have been created in an encounter be¬ 
tween two major galaxies in the distant past, a scenario 
which might explain the formation of thin planes of galax¬ 


ies such as the one observed around Andromeda I Pawlowski 
et al. 2011 Hammer et al.|2013 1 and around the Milky Way 


(e.g., Kroupa et al.|2005 2010 Pawlowski et al.|2012 |. Such 
tidal dwarf galaxies would be mostly devoid of dark mat¬ 


ter (Kroupa 19971, and would presumably be much more 


prone to tidal deformation than those embedded in primor¬ 
dial dark matter halos. Some of the problems of this picture 
are in principle alleviated if the Modified Newtonian Dynam¬ 


ics framework (Milgrom 1983 Famaey & McGaugh 2012 


Zhao et al.|2013 ) is adopted ( Kroupa et al.|2(5To Pawlowski 

& McGaugh|2014 1. Some other solutions consider the possi¬ 
bility of dissipative dark matter, as for example, in [Foot ^ 
Silagadze (20131 or in Randall & Scholtz (20141. 


In cosmological simulations, subhalos are generally 
rounder than the host halos. Subhalos that have masses that 
are large enough to host dwarf galaxies have axis ratios : 
0.65 ^ c/a ~ h/a ^ 0.95 ( Vera-Ciro et al.]|2014 (. Observed 
dSphs are observed to be approximately roughly round: for 
the Milky Way, e = 1 — a/b ~ 0.3 ( Martin et al.|[2008 (. But 
we know that this observed ellipticity is just a projected 
view, and the third dimension could thus hide more elon¬ 
gated shapes. However, deprojecting to recover the intrinsic 
ellipticity is a challenge because distance uncertainties to 
different sides of the same dSph are currently far too large 
to allow any differential depth measurement except maybe 


in the Magellanic Clouds (Moretti et al. 2014 Klein et al. 


2014). Therefore, in this paper, we propose a method which 


allows us to probe the third dimension, and we focus our¬ 
selves on Mof] 

In sectionj^we summarize the general method employed 
in this contribution, justify our hypotheses and present our 
data. Section presents all dwarf galaxies parameters used 
in the analysis, which have been extracted from the observa¬ 
tions. In section[^ we define how our deprojection procedure 
is built. Results are then discussed in sections and and 
our conclusions are drawn in section [3 


shape of dark matter halos in the relevant mass range (Vera- 
Ciro et al.|2014 Barber et al.|2015 |. In simulations in which 
dwarf satellites galaxies are initially made of rotating stellar 
discs, this is a less natural assumption, but they can in prin¬ 
ciple also be deformed into a prolate shape through bar-like 
instabilities caused by ongoing small tidal interactions with 
the host ( Lokas et al.|2010 |. Our second assumption will be 
to take an orientation of the prolate body that is consistent 
with the observed projection of the major axis, as sketched 
in Figure Three cases will be considered: 

(i) the prolate dwarf galaxy is oriented perpendicular to 
the radial vector linking the host and the satellite 

(ii) the major axis of the prolate dwarf galaxy is oriented 
as close as possible to the radial vector 

(iii) the prolate dwarf galaxy has a random orientation 
that is consistent with the observed projection. 


2.2 Data 

In the present contribution, we will examine the struc¬ 
tural properties of the dwarf galaxies surrounding the An¬ 
dromeda galaxy (M31). Our observational data are derived 
from the Pan-Andromeda Archaeological Survey (PAndAS) 


(McConnachie et al. 2009). This photometric survey with 


the MegaCam camera on the Canada-France-Hawaii Tele¬ 
scope has mapped a large area around Andromeda (~ 400 
deg^), allowing the discovery and characterization of the 
dwarf satellite galaxies of M31. A detailed account of the 
processing of the PAndAS survey through to the final cata¬ 
logue of aperture magnitudes is given in ( Ibata et al.|2014 1. 

The PAndAS photometry was also used to estimate the 
distances to the dwarf galaxies, using the tip of the red gi¬ 
ant branch (TRGB) method ( [Conn et al.|2012 ). By using a 
Markov Chain Monte Carlo (MCMC) approach that incor¬ 
porated Bayesian priors on the expected morphology of the 
satellites and their luminosity function, [Conn et al.| ( |2012[ ) 
were able to derive a posterior probability distribution func¬ 
tion (PDF) for the distance to each of the satellites present 
in the survey. 

From the homogeneous PAndAS data set, we also in¬ 
ferred the morphological parameters of all dwarf galaxies 
that fall in the survey footprint. These were derived using a 


version of the Martin et al. (2008) algorithm, updated for a 
full MCMC treatment, and which will be presented in an¬ 
other publication that is to be submitted soon (Martin et 
al. 2015, in preparation). 


2 METHOD 
2.1 Hypotheses 


3 EXTRACTING DWARF GALAXY 
GHARAGTERISTICS FROM THE 
OBSERVATIONS 


In order to deproject the view and try to see which intrinsic 
ellipticity distribution can fit the observed projected distri¬ 
bution of stellar populations in the target dwarf, some sim¬ 
plifying hypotheses need to be made. In the following, two 
main hypotheses are adopted. First, we consider that dSphs 


^ The method presented here could be modified later on to be 
applied to dwarf satellites of the Milky Way. 


3.1 Morphology and selection 


The code developed by Martin et al. (20081 infers the fol¬ 
lowing observational parameters : the projected ellipticity 
(e, with e = 1 — (a/b) where a is the minor axis of the ad¬ 
justed ellipse and b the semimajor axis), the position angle 
(PA, the angle between the equatorial North and the major 
axis), and the half-light radius (r^, taken along the major 
axis). The code was run on 23 dwarf galaxies present in the 
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Plane of the sky 


Plane ‘P defined by the line of sight (los) and 


Projected view 


North 




Figure 1. Sketch of the geometry. The left-hand panel shows the projected view on the sky, while the right-hand panel shows the view 
in the plane P defined by the line of sight and the projected semimajor axis. Two cases are represented : the tangential case and the 
radial one. The bottom right-hand panel shows the position of the plane "P compared to M31 where /3 is the minimum angle between the 
radial direction and P. The variable a is the semi-minor axis of the prolate ellipsoid, equal to the semi-minor axis of the projected ellipse. 
On the left-hand panel, PA is the angle between North and the semi-major axis ff^ where ||rj^|| is the observed half-light radius. On 
the right-hand panel, hr is the intrinsic semi-major axis of the prolate ellipsoid placed in the configuration that is most radially-oriented 
with respect to M31, and bt is the one placed in the tangential configuration. The variable 7 is the angle between the major axis of the 
prolate ellipsoid and its intrinsic major axis. The line of sight position vector los links the Earth to the centre of the satellite galaxy. 


PAndAS survey and outputs posterior PDFs in the form 
of MCMC chains. We have also included in our study the 
dwarf elliptical satellites, NGC 147 and NGC 185. For them, 
realistic MCMC chains have been created by simultaneous 
drawings from gaussian PDFs built from [Crnojevic et ah] 
(20141, also derived from the PAndAS survey. According to 


Ibata et al. (20131, among these 25 satellites, 14 are in the 


VTP: And I, And III, And IX, And XI, And XII, And XIII, 
And XIV, And XVI, And XVII, And XXV, And XXVI, 
Cass II, NGC 147 and NGC 185. Whereas 11 are outside 
the VTP: And II, And V, And X, And XV, And XVIII, 
And XIX, And XX, And XXI, And XXII, And XXIII and 
And XXIV. 


Thus, for each satellite k of Andromeda, the code pro¬ 
vides us with a Markov chain that contains at each step 
i the positional and morphological parameters of interest: 
Chi, PAm and For example, in Figure[^we display the 
PDF derived from 10® draws (this is typically the number 
of draws in all MCMC chains used in this study) of the 
projected ellipticity from the Markov chain of Andromeda 
I, an example of a rather round dSph galaxy with a large 
number of resolved red giant stars and Andromeda XXII, an 
example of an extended dSph galaxy with a small number 
of resolved red giant stars. It is more representative to plot 
log(l — e) = log(a/b) where 0 corresponds to a round object. 

When the whole sample of dSphs considered is taken 
into account, one draw from their respective Markov chain 


of ellipticity and distance represents a possible view from 
the Earth of the system of satellites around Andromeda. 
Figure [^represents the normalised average of the projected 
ellipticity of the satellites according to the derived PDFs 
for 10® random realisations (of the 25 satellites). Bins have 
been chosen to be large enough, log{l — e) — 0.2, to ensure 
reliable statistics. The vertical bars are one sigma standard 
deviations of the number of satellites per bin, derived from 
7 X 10"^ draws. 

As expected, the projected ellipticity shows predomi¬ 
nantly round objects. Moreover, no signihcant difference be¬ 
tween galaxies in the VTP (blue triangles) or outside it (grey 
squares) is evident. 


3.2 Position in three dimensions 


Once morphological parameters have been obtained, the po¬ 
sition of each dSph centre in three dimensions has to be 
found. For that purpose, Andromeda’s satellite galactic co¬ 
ordinates (ofc, Sk) are considered perfect. We have chosen 
to take them from the literature in Collins et al.j ( |2013 1 for 
Cassiopeia II and in |McConnachie| ( |2012[ ) for the other ones. 
For the third dimension, the distance from Earth (d), we use 
results found by Conn et al. (20121 (see Figure]^. This pro¬ 
vides us with a Markov chain of distance for each galaxy dk- 
Each step i of it gives a possible dki for the satellite. 

For a dwarf galaxy, a drawing in its morphological 
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Figure 2. Representative projected ellipticity distributions. His¬ 
togram of the projected ellipticity for And I (solid line) and And 
XXII (dotted line) (two representative dSph in the sample) ob¬ 
tained from 10^ draws from the Markov chain of Martin et al. 
(2015) (in preparation). 


. Galaxies inside the VTP ^2d') 



log(l-e) 


Figure 3. Observed (i.e., projected) ellipticity derived from the 
PAndAS survey, showing the average distribution of ellipticity for 
the 25 satellites, randomly drawn 10^ times from each of the 25 
Markov chains. Bin sizes are log{l — e) = 0.2. Red circles are the 
ellipticity distribution of the entire system of 25 dwarf galaxies. 
Blue triangles correspond to the 14 dwarf galaxies contained in 
the vast thin plane of satellites ( |Ibata et al.|201^ , while the grey 
squares are for the 11 outside of the VTP. Vertical error bars 
indicate one standard deviation derived from this experiment. 


Markov chain and one in its distance Markov chain com¬ 
bined with sky coordinates provide a set of six parameters. 
They define a possible representation in space of the centre 
around Andromeda in association with three possible pro¬ 
jected morphological characteristics : Sk, dki, €ki, PAki 

and 


Figure 4. Histogram of the distance of the centre of two repre¬ 
sentative systems, And I (solid line) and And XXH (dotted line) 
from the centre of M31 (derived from 10^ draws from the Markov 
chain of |Conn et al. |2012^ . 


3.3 Model 


In all the following analysis, the intrinsic shape of the dwarf 
galaxies is approximated as a prolate ellipsoid, which is a 
natural outcome of dark matter satellites in ACDM simula¬ 
tions i Vera-Ciro et al.]|2014 Barber et al.|[2015 l. With this 
assumption the satellites have two equal minor axes (a) and 
one major axis (6). Consequently, since a prolate structure 
is symmetrical by rotation around the major axis, the pro¬ 
jected minor axis remains identical, irrespective of the point 
of view. Thus the apparent minor axis of the ellipse given 
by the projected satellite is taken to measure the real minor 
axis (a). On the other hand, fh is known, defined by PA 
and Vh- 


Of course, due to the projection on the sky, the real 
major axis b is not equal to rTl: b can have any orientation 
in the plane {P) defined by the line of sight and rjj (see Fig¬ 
ure [^. Thus, an additional condition is required to fix the 
three-dimensional orientation of the major axis, and conse¬ 
quently the intrinsic ellipticity of the satellites. We consider 
three possibilities which will be discussed in Section]^ 


3.4 Limitation of the major axis 

As defined above, the intrinsic major axis b has to lie in the 
plane P. However, depending on the parameter set drawn 
from the chain, it is possible that b can take on physically 
implausible values. To avoid this situation, we implement a 
filter to reject solutions where b would be larger than the 
tidal radius of the satellite (tidal rejection criterion). Thus, 
the tidal radius is considered as the maximal radius : 



where Msat is the satellite mass and Mhost the mass of the 
host galaxy contained in the sphere of radius R (distance 
between the satellite centre and the host galaxy centre). We 
estimate Mhost using an NFW (Navarro et al. 19971 halo 
model (with virial mass 2 x 10^^ Mq, a mean value of recent 
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results ( Watkins et al.|2010 Fardal et al. [2013 1). The value 
of Msat is derived from the line of sight velocity dispersion 
((j„) listed in McConnachie (|2012 I for NGC 147 and NGC 
185 and in Gollins et al. ( 2013[) forthe other galaxies, taking: 


Ms, 


( 2 ) 


where G is the gravitational constant. Gonsequently, with 
equation 0 and 0 . the following limit can be set to b : 


b < (Jv 




2 G Mhost 


( 3 ) 


Given the large uncertainties in (t„, we wish to remain con¬ 
servative and so as not to eliminate physically possible cases, 
we adopt the values corresponding to the listed value of cr„ 
plus one standard deviation. In the following analysis, this 
criterion is always automatically applied, and it forces the 
rejection of approximately 5% of the MCMC parameter trial 
values. Finally, to avoid the MGMG process from exploring 
unrealistic solutions, we limited the projected half-light ra¬ 
dius to 20 arcmin. 


3.5 Simple test 

We constructed a simple toy model to test the algorithms 
that we developed and to obtain a first assessment of the 
magnitude of the effect of de-projection on the ellipticity 
distribution. Three-dimensional prolate dwarf galaxies are 
generated with a random size (in the range of the observed 
dwarf satellites galaxies) and are placed at random locations 
around the Andromeda galaxy (again in the observed range). 
In Figure we show the consequence of assuming that the 
intrinsic ellipticity is given by the observed projected distri¬ 
bution previously displayed in Figure (red dots). These 
prolate structures are observed from a position correspond¬ 
ing to Earth, and we measure the projected ellipticity that 
the artificial structures would then be seen to possess (purple 
triangles). As expected, the projected distribution will ap¬ 
pear rounder than the intrinsic distribution. The differences 
in the rounder bin is about 2a, thus confirming that the 
projected distribution cannot be seen as the 3-dimensional 
one. 


Projection test 



- 0.8 - 0.6 - 0.4 

log(l-e) and log(l—E) 


Figure 5. Green circles represent the intrinsic ellipticity cho¬ 
sen for this test. Purple triangles represent the distribution of 
projected ellipticity for 1000 3-dimensional systems of 25 dSph 
galaxies generated from a random distribution of position and 
from the green distribution for the ellipticity. Red dots are the 
observed ellipticities previously seen in Figure Vertical error 
bars indicate one standard deviation. 
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4 DEPROJECTION 

Our next task is to hnd the distribution of intrinsic ellip¬ 
ticities that is consistent with the observed projected distri¬ 
bution of Figure]^ and with the observed orientations and 
distances of the real satellites. As stated above, this will 
require an additional geometric constraint, which we will 
implement via an assumption regarding the 3-dimensional 
orientation of the satellites. 

4.1 Tangential deprojection 

As a first hypothesis, we suppose that each satellite is ori¬ 
ented tangentially with respect to the vector fj, connecting 
it to M31, i.e. that the real semimajor axis (6) of the mod¬ 
eled dwarf galaxy is on the plane V' perpendicular to ft,- 
For satellite k, a given set i of parameters {Ik, bk, dki) then 


Figure 6. Histogram of the ellipticity for And I (solid black line) 
and And XXII (dotted black line) after tangential deprojection 
obtained from 10® draws of the MCMC chain. The red curves are 
the observed ellipticity previously seen in Figure]^ 

uniquely defines the direction of the vector b, which corre¬ 
sponds to the intersection of V and V' ■ Once the direction of 
b is known, 116| | is determined by requiring that it correspond 
to the projected view of the dwarf galaxy whose apparent 
semimajor axis is ||r7,||. Finally, the intrinsic ellipticity Et^. 
assuming this “tangential hypothesis”, is calculated (see Fig¬ 
ure]^. In Appendix [A| we describe this deprojection in more 
detail. 

As the tangential deprojected chain of each of the dwarf 
galaxies is now known, a random draw from each of the 
Markov chains gives a system of satellites all of which are 
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Perfectly tangential 

0.8 ^ 0.5 0.2 



Figure 7. Distribution of tangentially deprojected ellipticity for 
10® systems generated from Markov chains derived from the ob¬ 
servational analysis. Green circles mark the ellipticity distribution 
of the entire system of 25 dwarf galaxies. Blue triangles are for 
those contained in the vast thin plane of galaxies pbata et al.| 
|2013| l, 14 dSphs and grey squares are for those outside the VTP, 
11 dSphs. Red dots are the observed ellipticities previously seen 
in Figure]^ Vertical error bars indicate one standard deviation. 

in tangential orientation. This procedure is repeated to ob¬ 
tain the global ellipticity distribution of the dwarf galaxy 
system, which is shown in Figure Under this “tangential” 
assumption, the intrinsic distribution of ellipticity is sub¬ 
stantially flatter than the projected ellipticity, and produces 
a long tail of highly flattened systems. No significant dif¬ 
ference between galaxies within the VTP or outside of it is 
apparent. 

4.2 Radial deprojection 

The second hypothesis that we explore in this analysis is 
to consider that the satellites are oriented as close as pos¬ 
sible to the radial direction with respect to Andromeda. As 
before, the major axis b must lie on V. For this “radial” 
deprojection, we would like b to lie along rt,- However, the 
prolate ellipsoid must also be consistent with the observed 
(projected) orientation of the galaxy. In this situation, we 
have three spatial unknowns that have to respect four inde¬ 
pendent equations, so that the problem turns out to be over¬ 
constrained. Thus, an exact radial orientation is, in general, 
impossible. 

A way to resolve the problem is to attempt to orient the 
prolate ellipsoid as close as possible to a radial direction. 
To this end, we minimise the angle /3, between and b 
such that the projection of b yields rj). Consequently, j3 is 
the minimum angle, an incompressible lower limit, that can 
exist between the radial vector and the major axis vector of a 
dwarf galaxy according to our data. The intrinsic ellipticity 
(assuming this “radial” hypothesis) is thus obtained 

(see Figure]^ coupled with /3m (see Figure]^. For more 
details, see appendix [B| 

The radially deprojected ellipticity is displayed in Fig¬ 
ure]^ and interestingly, the roundest bin is not dominant. 


0.8 0.5 0.2 



Figure 8. Histogram of the ellipticity for And I (solid black 
line) and And XXII (dotted black line) after radial deprojection 
obtained from 10® draws of the MCMC chain. The red curves are 
the observed ellipticity previously seen in Figure]^ 



Figure 9. Probability density function of the angle /3 between 
the major axis of : And I (solid line), And XXII (dotted line) 
and their respective radial direction after radial deprojection (ob¬ 
tained from 10® MCMC draws). 

and a significant number of highly elliptical systems are re¬ 
quired if this assumption is realistic. Again, differences be¬ 
tween galaxies within the VTP of galaxies and outside of it 
are not significant. 

4.3 7 -deprojection (random angle within V) 

In the previous subsections we considered a “tangential” and 
a “radial” de-projection. We have seen that a satellite can 
always be fit to a 3-D tangential orientation, but that a 
prolate system cannot be placed closer to a radial direction 
than a minimum angle (3 (which depends on the structure 
and position of the satellite). However in reality, we do not 
have observational information that allows us to favor ei- 
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Figure 10. As Figure]^ but for the radial deprojection method. 


Figure 12. Probability density function of the angle 7 (in de¬ 
grees) between the major axis of And I (solid line) and And 
XXII (dotted line) with their respective radial direction after 7 - 
deprojection (obtained from 10^ MCMC draws). 
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Figure 11. Histogram of the ellipticity for And I (solid black 
line) and And XXII (dotted black line) after 7 -deprojection ob¬ 
tained from 10^ draws of the MCMC chain. The red curves are 
the observed ellipticity previously seen in Figure]^ 


Figure 13. As Figure but for the 7 deprojection method 
(where we assume that the major axis lies in plane 7^ at a random 
angle with respect to the line of sight). 


ther of these orientations, and we only know that the major 
axis lies within the line of sight plane V. This motivates our 
final deprojection, where we instead make the minimal as¬ 
sumption that the satellites are randomly oriented within V. 
The random angle 7 , is the angle between the line of sight 
and the major axis b. The intrinsic ellipticity (assuming this 
“gamma” hypothesis) ^~tki is thus obtained (see Figure [Tl| 
coupled with 7 ^; (see Figure [7^. 

The resulting deprojection is shown in Figure and 
presents a rather flat ellipticity distribution. Again, there is 
almost no difference between the dwarf galaxies that belong 
to the VTP and those outside of it. 


5 RE-PROJECTION TESTS 

The distributions of ellipticity obtained in Figures [7| and p0| 
take into account the observational information with a hy¬ 
pothesis on the orientation and on the prolate shape of the 
satellite galaxies. In contrast, the distributions of ellipticity 
obtained in Figure p^ take into account the observational in¬ 
formation assuming that the shape of the satellite galaxies 
is prolate, but that no particular direction can be privileged. 
If either of these approaches is realistic, one of the resulting 
deprojected distributions should be the real intrinsic ellip¬ 
ticity. 

To test the consistency of this approach, we decided to 
examine whether the derived deprojected ellipticity distri- 
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Perfectly tangential 

0.8 ^ 0.5 0.2 



log(1—E^) , log(1—e^) and log(l-e) 


y oriented in plane P 

0.8 ^ 0.5 0.2 



Figure 14. Distribution of projected ellipticity for 1000 sets of 
25 dSph satellites randomly built with intrinsic ellipticity picked 
from the tangentially deprojected distribution (purple triangles). 
Red dots show the distribution of observed ellipticities. Green 
circles are the intrinsic ellipticity determined with the tangential 
deprojection method represented in Figure]^ Vertical error bars 
indicate one standard deviation. 


As close to radial as possible 

0.8 ® 0.5 0.2 



Figure 15. As Figure [Tl] but for the radial deprojection method. 


butions would give rise to projected distributions similar to 
the observed distribution of e for random viewing directions. 

To this end we modelled systems containing 25 (3- 
dimensional) prolate ellipsoids with random orientation. The 
ellipticity, and so the semimajor axis, is randomly picked 
from the 3-D tangential-deprojected ellipticity distribution. 
The positions of the dwarf galaxy models are drawn ran¬ 
domly inside a sphere of 450 kpc around the host. A set of 
1000 such systems were generated and observed from a po¬ 
sition equivalent to that of the Earth point at 780 kpc from 
the host galaxy. The final 2-D projected ellipticity is shown 
Figure |14| The same test is applied with the 3-D radial- 
deprojected ellipticity and the 3-D 7 -deprojected one. Re¬ 
sults are shown Figures and respectively. 


Figure 16. As Figure [I^ but for the 7 -deprojection method. 


In Figure one can see that the projected ellipticity 
distribution (purple triangles) based on the radial deprojec¬ 
tion does not reproduce well the observed distribution of 
ellipticity (red circles) in the roundest bin. This means that 
if dwarf galaxies are prolate structures, and if our view¬ 
ing direction is not special, the galaxies are not oriented 
in a radial conhguration. However, the projected distribu¬ 
tions based on tangential (Figure [T^ or the 7 -deprojection 
(Figure |16[ ) methods are perfectly consistent with the ob¬ 
served distribution of ellipticity. We consider that this is a 
good consistency test of our method. Thus modelling dwarf 
galaxies as possessing prolate shapes can give self-consistent 
solutions in terms of their derived ellipticity distribution. 
Furthermore, the population of dwarf galaxies around An¬ 
dromeda can be considered as a drawing from a distribution 
of general intrinsic ellipticity. 


6 INDIVIDUAL CASES AND DISCUSSION 


The results for all of the dSph galaxies are listed in table 
The second column contains the character ”p” if the galaxy 


belongs to the VTP according Ibata et al. (20131, and ”a' 


otherwise. Columns 3 and 4 contain the mean observed el¬ 
lipticity and the mean observed half-light radius obtained 


with the method of Martin et al. (20081. Column 5 is the 


mean distance from M31 to the centre of the dSph galaxy 
obtained by jConn et al. (2012 1 . Note that these observed av¬ 
erage values (columns 3, 4 and 5) take into account the fact 
that some values have been rejected by our tidal criterion in 
the 7 -deprojection method. Columns 6 , 8 and 10 represent 
respectively the mean of the tangential-deprojected, radial- 
deprojected and 7 -deprojected ellipticity obtained with our 
method. Columns 7, 9 and 11 are the percentage of excluded 
cases by our tidal criteria for the 3 deprojections. The last 
column is the angle /? previously defined as the geometri¬ 
cally minimum possible angle between the radial direction 
and the plane V. All uncertainties indicate one standard 
deviation of each parameter with our method. 

First, the most powerful geometric constraint is a high 
angle /?. Indeed, it corresponds to the lower limit for a ra- 
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Table 1. Properties of the 25 sate llites of Andromed a studied here, s howing the results of th e 3 different deprojections obtained from 
10® draws of the MCMC chain. (l) |lbata et al.| ( |2013^ (2) According t o|Crnojevic et al.| l |2014| | for NGC147 and NGC185 and to Martin 
et al. (2015) (in preparation) for the other galaxies (3) [Conn et al.| ( |2012^ * Note that these values are calculated according to the 
7 -deprojection tidal rejection criterion. 


Name 

VTP 

(1) 

e 

(2)* 

(arcmin) 

(2)* 

dhiil 

(kpc) 

(3)* 

Et 

Tide 

(%) 

Er 

Tide 

(%) 


Tide 

(%) 

P 

(deg) 

Audi 

P 

0.29±0.03 

3.97±0.14 

68.7±13.3 

0.37±0.05 

1.8 

0.51±0.10 

0.0 

0.43±0.15 

14.3 

26.7±6.0 

Andll 

a 

0.14±0.02 

5.13±0.10 

197.0±11.0 

0.15±0.02 

0.0 

0.57±0.09 

0.3 

0.31±0.20 

5.0 

43.4±3.3 

Andlll 

p 

0.59±0.03 

1.89±0.15 

91.2±13.9 

0.66±0.05 

0.4 

0.76±0.06 

0.0 

0.71±0.13 

6.8 

40.9±7.7 

AndV 

a 

0.28±0.07 

1.68±0.14 

116.3±5.9 

0.70±0.15 

7.9 

0.31±0.07 

0.0 

0.47±0.22 

3.3 

16.2±7.1 

AndIX 

p 

0.12±0.07 

1.85±0.21 

138.4±52.3 

0.12±0.08 

0.0 

0.63±0.23 

6.5 

0.29±0.23 

1.9 

11.0±7.2 

AndX 

a 

0.30±0.18 

1.15±0.24 

143.9±27.6 

0.36±0.20 

0.1 

0.54±0.19 

0.4 

0.47±0.28 

2.6 

12.5±8.9 

AndXI 

p 

0.24±0.16 

0.78±0.20 

141.0±41.7 

0.46±0.31 

2.6 

0.41±0.27 

0.4 

0.43±0.29 

1.3 

34.4±24.1 

AndXII 

p 

0.39±0.25 

2.85±1.36 

157.7±52.2 

0.47±0.26 

15.2 

0.56±0.24 

5.1 

0.50±0.27 

18.7 

15.9±12.9 

AndXIII 

p 

0.54±0.16 

1.11±0.41 

180.2±81.2 

0.71±0.19 

6.9 

0.65±0.17 

0.0 

0.68±0.19 

3.1 

9.8±7.4 

AndXIV 

p 

0 .20±0.11 

1.55±0.16 

184.4±34.9 

0.58±0.28 

10.1 

0.25±0.13 

0.0 

0.39±0.25 

3.0 

13.4±10.4 

AndXV 

a 

0.24±0.10 

1.39±0.13 

167.9±45.5 

0.26±0.10 

0.4 

0.73±0.20 

18.7 

0.42±0.24 

3.5 

36.2±13.7 

AndXVI 

p 

0.29±0.08 

l.OOiO.08 

316.1±33.7 

0.30±0.08 

0.0 

0.79±0.08 

0.2 

0.49±0.23 

0.7 

21.9±3.0 

AndXVII 

p 

0.42±0.10 

1.39±0.21 

71.8±21.0 

0.49±0.12 

7.4 

0.66±0.15 

0.5 

0.55±0.16 

13.8 

30.5±11.0 

AndXVIII 

a 

0.21±0.13 

0.78±0.11 

454.6±40.1 

0.22±0.14 

0.0 

0.73±0.14 

0.4 

0.40±0.28 

0.4 

7.6±4.8 

AndXIX 

a 

0.41±0.09 

11.65±1.64 

183.8±55.2 

0.49±0.08 

74.4 

0.56±0.17 

65.2 

0.49±0.12 

72.5 

3.7±2.7 

AndXX 

a 

0.30±0.20 

0.77±0.25 

141.7±13.4 

0.57±0.31 

2.9 

0.38±0.24 

0.2 

0.48±0.30 

1.5 

32.5±24.9 

AndXXI 

a 

0.32±0.11 

4.18±0.61 

137.2±9.7 

0.54±0.14 

12.9 

0.39±0.12 

0.1 

0.44±0.17 

16.9 

36.2±10.8 

AndXXII 

a 

0.58±0.14 

1.14±0.32 

261.5±33.6 

0.77±0.13 

3.3 

0.63±0.13 

0.0 

0.72±0.17 

3.2 

17.1±8.6 

AndXXIII 

a 

0.38±0.05 

5.43±0.40 

130.3±3.4 

0.72±0.07 

57.4 

0.40±0.06 

0.0 

0.51±0.14 

16.1 

35.3±4.7 

AndXXIV 

a 

0.23±0.15 

2.92±0.72 

166.7±24.8 

0.30±0.19 

0.9 

0.45±0.24 

2.9 

0.38±0.26 

6.8 

25.8±13.4 

AndXXV 

p 

0 .20±0.11 

3.18±0.33 

111.9±32.2 

0.33±0.17 

21.2 

0.32±0.17 

1.3 

0.31±0.18 

18.3 

23.3±14.0 

AndXXVI 

p 

0.30±0.20 

1.45±0.49 

197.6±108.9 

0.43±0.27 

2.7 

0.54±0.28 

0.8 

0.48±0.29 

2.1 

20.6±17.0 

CasII 

p 

0.38±0.07 

1.47±0.14 

176.2±40.4 

0.45±0.09 

0.3 

0.69±0.12 

0.0 

0.57±0.21 

1.1 

36.1±10.0 

NGC147 

p 

0.46±0.02 

6.70±0.09 

118.6±10.1 

0.63±0.07 

1.6 

0.58±0.06 

0.0 

0.60±0.15 

10.4 

39.1±4.8 

NGC185 

p 

0 .22±0.01 

2.94±0.04 

181.9±15.7 

0.24±0.01 

0.0 

0.60±0.04 

0.0 

0.42±0.23 

1.2 

24.8±2.3 


dial orientation. This means that even if the orientation is 
unknown, a dwarf galaxy cannot be closer to the radial di¬ 
rection than P (see Figure 1). This criterion does not depend 
at all on our hypothesis. Thus, in the present study, we have 
identified 6 dwarf galaxies - And II, And III, And XXI, And 
XXIII, CasII and NGC 147, which can absolutely not be ra¬ 
dially oriented. For this we have chosen a criterion such that 
P — 1g > 25° (where cr is here the standard deviation of P). 
These dwarf galaxies are relatively close to their host. More¬ 
over, they mostly have a large half-light radius and favour 
a high ellipticity (more than 0.5) in each kind of deprojec¬ 
tion. Also, for 5 of these 6 galaxies, the 7 -deprojection shows 
the largest rate of tidal rejections among the 3 cases. These 
particularities suggest that these six dwarf galaxies may be 
largely affected by their environment and could be affected 
by tides or even remnants of ancient tidal dwarf galaxies 
(TDG). For the most massive object, NGC 147, This last 
result is in a good agreement with Crnojevic et al. (20141 
where they have recently found that NGC 147 presents tidal 
tails. Moreover, it could also indicate that it is beginning to 
be stretched due to gravitational interactions, as attested 


by the stellar stream that emanates from it (McConnachie 
|et al.|2009| ). 

It is also interesting to see the deprojection behaviour 
of the other dwarf galaxies that have an observed large half- 
light radius and a high level of tidal rejections with all 3 
deprojection methods. And XII, And XIX and And XXV 
are such cases. This might be due to the effects of tidal in¬ 
teractions with M31, as argued by Collins et al. (20141 for 


And XIX and And XXV (as well as And XXI). Moreover, 
this strongly suggests that these galaxies are close to being 
destroyed and becoming stellar streams. Please note that 
we have been conservative in estimating the tidal limits of 
the dwarf galaxies, as we used the measured velocity dis¬ 
persion plus one standard deviation to calculate the tidal 
limits. Thus the percentage of tidal rejections listed in Ta¬ 
ble [T] should be considered as lower limits. 

Finally, And XIII and And XXII have a large intrin¬ 
sic ellipticity with all of the deprojection choices. Indeed, 
with the 3 deprojections, they always conserve an ellipticity 
around 0.7. Thus, these dSph galaxies are intrinsically very 
elongated. 

As the re-projection tests suggest a possible tangen¬ 
tial orientation being very plausible, we have examined in 
this case the deprojected ellipticity as a function of the dis¬ 
tance (Figure [T7| ). Blue circles are the values for the satellites 
within the VTP and grey circles are the values for those out¬ 
side. There is a priori no obvious correlation. Nevertheless, 
a close inspection of the repartition of the satellites outside 
of the VTP seems to indicate a relatively tight correlation 
between 100 and 200 kpc : dSphs are more elongated when 
closer to M31. No such correlation is present for satellites in 
the VTP. With the large error bars, no strong conclusion can 
be drawn, but it would be interesting to compare this with 
different models, both in standard cosmology in which the 
VTP could have been recently accreted and in alternative 
scenarios such as those based on TDGs. 

We next examine the properties of the population of 
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Figure 17. Tangentially deprojected ellipticity versus distance 
to Andromeda. Blue circles are the values for the satellites in 
the VTP of galaxies pbata et al.|201^ while grey circles are for 
the ones outside. The filling of these circles represents the tidal 
rejection rate (see labels in the top-right box). 


Table 2. Average ellipticities of the whole sample of dwarf galax¬ 
ies of Andromeda used in this work, obtained from 25 X 10® draws 
of the MCMC chain. 



Mean ellipticity 

Tidally rejected 
(%) 

e 

0.32±0.11 

— 

Et 

0.44±0.14 

9.22 

Er 

0.54±0.15 

4.12 

E-y 

0.48±0.21 

9.14 


25 dwarf galaxies as a whole, which we consider as a rep¬ 
resentative sample for this kind of object. The ellipticity of 
this group is listed in Table for projected (e), tangential- 
deprojected (Et), radial-deprojected (Er) and 7 -deprojected 
(E.^) with their total percentage of tidal rejection. This per¬ 
centage is compatible with a physical behaviour with a high 
value for tangential and 7 deprojections and a much lower 
value for the radial deprojection. Objects in a radial config¬ 
uration are falling toward their host with a compact shape 
whereas objets in tangential or random orientation are sus¬ 
pected to have been strongly disrupted in their orbits around 
their host. 

As expected, each kind of deprojection gives a larger 
ellipticity (with a minimum of 0.44) than the projected ob¬ 
servations. Thus, these dwarf galaxies tend to possess an 
elongated shape with a minor to major axis ratio (a/b) of 


approximately 1/2. Vera-Ciro et al. (20141 and Barber et al. 


(20151 find a sightly larger but coherent axis ratio in the 
inner regions (~ 1 kpc) of the dark matter satellites in the 
Aquarius simulation that could possibly host stars. Thus our 
conclusion appears to be consistent with the possibility that 
stellar components in the Andromeda satellite galaxies pos¬ 
sess the same shape. This is perhaps not incompatible either 
with them being remnants of ancient tidal dwarf galaxies. 


These observational values can in any case serve as a test of 
various models in the future. 


7 CONCLUSIONS 


First of all, we have seen that adopting a prolate shape for 
the dwarf satellite galaxies of Andromeda gives a coherent 
solution in the sense that the derived ellipticity distribution 
for the real viewing direction is similar to what would be 
derived from a random vantage point. This simplifying as¬ 
sumption has allowed us to derive the intrinsic ellipticity for 
the entire population of satellites. The resulting ellipticity 
distribution is found to be relatively flat between log{l — e) of 
0 to —0.4 before decaying at larger ellipticities. Some of the 
dwarf galaxies are found to be intrinsically significantly elon¬ 
gated. This is probably not incompatible with them being 
remnants of ancient tidal dwarf galaxies as recently found 


by, e.g., Yang et al. (20141, but it could also be consistent 


with dwarf galaxies having the same shape as their hosting 
dark matter subhalo. Secondly, no clear difference in elliptic¬ 
ity between the satellites inside the VTP and those outside 
of it was noticed. The two populations of galaxies seem to 
be made of the same type of objects as also recently con¬ 


cluded by Collins et al. (20151. Nevertheless we noted that 


the distribution of ellipticity versus distance to M31 in the 
tangential deprojection case was slightly different for galax¬ 
ies inside and outside of the VTP. 
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APPENDIX A: TANGENTIAL CASE 


We seek the vector m which has the same projection on the 
sky as the half-light radius rjj. This means it has the same 
projected position angle (PA) (see eq. A31. At the same time, 
rfi must be perpendicular to the direction ft,, between M31 
and the centre of the dSph (see eq. |A2[ ). First of all, for a 
dSph, we define a local coordinate system {los,n,e) centred 


on the galaxy where the line of sight los is one of the axes 
and where n points towards the north. From (standard) sky 
coordinates (X:^) with respect to the centre, we obtain the 
local base: 


( cos ?7 sm X \ 

sin ?7 , ir = 

cos ?7 cos x/ 


- sm X sm rf \ 

cos rj 1 6 = 

- cos X sin 77 / 


cosx \ 

-smx/ 

(Al) 


From this base with our conditions and the normalisation of 
m, a 3 equation system can be deduced: 


rl.m = 0 


(A2) 


^ ^ m.e 

tan(P74) = ^ ^ 
m.n 


(A3) 


mj' + rriy^ + = 1 . 


The solution of this system gives us: 


-1 

B 
A 

\ VC 


i^+rsj^ 


where: 


(A4) 


(A5) 


A — tan(PA) cos x sin riVs,. — tan(PA)rs^ sin x sin rj 

—Ts^ COSX — J'sj, sinx 

B — tan(PA)rsj^ sin x sin rj + tan(PA) cos rj 
+rsy cos X 

(Avs 




# +1 


As we know fh, the norm of the tangentially deprojected 
semimajor axis bt has yet to be discovered. DSph are consid¬ 


ered as prolate (see eq. A 6 I. Thus, we are putting ourselves 
this time in the plane defined by los and m. The new local 
base (y,«) corresponds to the minor (a) and the semimajor 
axis (bt) of the real prolate. 


2 2 

^ + — = 1 
^ bT 


(A 6 ) 


For the present purposes, dSph around M31 are far enough 
to be considered as being at infinity. Gonsequently, Vf, is 
approximately seen between two parallel rays, los and los' 
with a separation of ur. We called a the angle between y and 
los. yo is the point at the intersection between y and los'. It 
follows that: 


sina = — 

yo 


From eq. |A 6 [ we obtain y and: 

djl _ 

dz 


V 


1 - 


(A7) 


(A 8 ) 


With eq. A 8 we obtain the equation for the straight line (d) 
which has los' for guiding vector: 


y = yo- 


1 - 


(A9) 


And so, the intersection between (d) (eq. A91 and the prolate 


(eq. A 61 with eq. A7 gives: 


1 - 


(AlO) 


The development of eq. with tana = —y' and /3t = 
^ — a gives the final real semimajor axis, tangentially de¬ 
projected: 


bt = 


sin Pt 


tan/ft 


(All) 
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APPENDIX B: RADIAL CASE 

We seek the vector q which has the same projection on the 
sky as the half-light radius rlj. It means the same projected 
position angle (PA). At the same time, q must be as close as 
possible to the radial direction fl, connecting M31 with the 
centre of the dSph (see eq. |A2[ |. We place the system in the 
same basis as the tangential deprojection (see appendix [A|). 
The line of sight is again los and m is the tangential direc¬ 
tion. Thus, q must be placed in the plane V defined by these 
two vectors. V has as its normal vector rPp: 

( 


n-p 


Bnp^ 

Anp^ 

1 


(Bl) 


\-/l+A^+SS , 


where: 


A = 
B = 




And so, with the minimum angle between q and ft,, the 
constraint of lying in Plane V and the normalisation of q, a 
system of 3 equations can be deduced: 


q.rfp = 0 

q.fs = max 

, 9^ + =1 


This system is solved by maximising eq. |B2| 

/'Fr(^r’,,qy+nv„qz)\ 

^ A --- 


9 = 


V 


‘ X 




where: 


(B2) 


(B3) 


= - nvyUvA - 

= - rivers,. 


Finding the norm br follows the same approach as that for 
the tangential case and so: 


br = 


( ^ 


\ sin pr 


tan P, 


(B4) 
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